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Abstract—In this paper, we propose a method based on signal 
processing techniques to improve the performance of aperture 
antennas, such as a parabolic reflector, used to estimate the 
Direction of Arrival (DOA) of a signal. Three signal processing 
algorithms are investigated, including the correlation method, 
which is used to make an initial estimate of the incidence angle 
within a certain range. Then inverse matrix method and singular 
value decomposition method are subsequently utilized to refine 
the estimate within this range. 


I. INTRODUCTION 


DOA estimation techniques for aperture antenna have wide 
applications in many fields, such as communication, radar, 
sonar, radio, etc [1]. Many algorithms, including the Multiple 
Signal Classification (MUSIC), Estimation of Signal 
Parameters via Rotation Invariance Techniques (ESPRIT), 
have been investigated to address this problem. However, these 
techniques do not perform too well on highly correlated or 
coherent signals generated via multipath propagation scenario. 
Alternative approaches for coherent signals have also been 
developed for DOA estimation, including the Maximum 
Likelihood (ML) method, Weighted Subspace Fitting (WSF) 
method, Spatial Smoothing (SS) method, Singular Value 
Decomposition (SVD) method and Toeplitz method, to name a 
few. Its accuracy of the ML method degrades severely when 
the angular separation between sources is small and the signal- 
to-noise ratio (SNR) drops below a certain level [2]. The WSF 
method presents a relatively high computational burden owing 
to the multi-variable minimization that it must perform. The SS 
method cannot resolve the number of signals that exceed the 
number of array sensors. The performance of SVD method is 
poor when the SNR level is low, because it reconstructs the 
matrix data set, which includes all the signal information, by 
directly processing the noise-free snapshot data[3]. In this 
paper, a novel technique for DOA estimation is proposed for 
aperture antennas, such as the parabolic reflector to improve 
the estimation in a computationally efficient manner. The 
proposed method uses the correlation approach to estimate the 
DOA within a certain angular range. The inverse matrix and 
SVD methods are subsequently employed to improve the initial 
estimate and improve the resolution of the aperture antenna. 
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II. | PROBLEM FORMULATION 


Fig.l shows a representative aperture antenna, namely a 
parabolic reflector. We assume that the sensors receiving the 
signals are distributed in the focal plane of the reflector along 
the y axis. Our task is to estimate the azimuthal angle of arrival 
by using the proposed method. 
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Fig. 1. Geometry of the parabolic reflector with sensors 
located in its focal plane. 


Uniform Linear Array (ULA) of N sensors receiving I 
signals from the directions 6;, i=1,*+,/. The incident field is 
assumed to be a plane wave originating from a large distance 
away from the antenna. The received signal vector can be as 
follows [4]: 


x(t) = As(t)+n(t) 
x() = [x (A), xy O] 
A =[a(0,),::,a(0,)] (1) 
s(t) = [s,s O 
n(t) = [7n 0, n (OV 
here x,(¢) represents the signal received at the nth sensor, with 
n=1,°*:, N; the matrix A comprises the vector of the signal 
waveform; s;(¢) is the ith signal, where i=1,-*-, I, and 7,(7) is the 
additive white Gaussian noise induced in the nth sensor. 
IHI. PROPOSED METHOD FORDOA ESTIMATION 


We begin the procedure, proposed herein, by using the 
correlation method as described below. The correlation 
coefficient p between two signals s;(f) and s,(f) is given by: 
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E[s, (1)s" («| 
Pik = 
Je [OF lello] 


where p represents the linear correlation (or dependence) 
between the two signals. If p=1, there is a total positive 
correlation while there is no correlation if p=0. 


(2) 


Following the correlation step, which provides as with as 
crude estimate of the range of the DOA, our next step is to 
investigate the SVD method, which is a type of matrix 
reconstruction technique that is useful dealing with ill- 
conditional matrices, which the matrix may well be for 
practical situations. For the noise-free case or if the SNR is low, 
we can use the Inverse Matrix method to find the incident angle 
of incoming signal. We note also that many DOA algorithms 
perform a subspace decomposition, by using an eigenvalue- 
eigenvector analysis, which is usually implemented by using 
the SVD method [5]. Rj; which is the output of the SVD 
operator, can be written as 


Rig = Cii Kapa (3) 


The normalized values show the positions of the incident 


angles, J is the number of the entries of R, C;,, is the inverse 


matrix of the new base matrix, which is the truncated version of 
the original base matrix, whose size N has been truncated to M 
by using the SVD, which delete the eigenvectors whose 
eigenvalues fall below a certain threshold, to improve the 
condition number of the new base matrix. The original base 
matrix can be derived in advance by using J different incident 
angles. We add that the Ky.,is the received signal at the NV 
sensors following the SVD projection. 


IV. NUMERICAL RESULTS 


For the purpose of the numerical experiments reported 
herein, we set the frequency of operation to be 30GHz, and use 
a reflector whose focal length (F) is 16.16mm, and it has an 
aperture diameter (D) of 63.6mm (6.364). We begin by using 
the HFSS simulator to generate the original base matrix for 61 
incident angles ranging from —30° to 30°, with a 1° interval. 
We assume that there are 67 sensors (0.0954, separation) on 
the focal plane along the y axis; hence N is 67. The SVD 


threshold is set to 10 ~, and the SNR is assumed to be 40dB. 


Fig.2(a) shows the magnitude distribution of the focal plane 
field when two incident waves from 0° and 3° are 
simultaneously incident upon the reflector. From the plot of the 
correlation coefficient shown in Fig.2(b), we see that we can 
narrow the search range within —2° and +5°. The DOA results 
for this reduced range are shown in Fig.3(b), and are compared 
with those in Fig.3(a) when the entire angular range is used. 
We can see that the proposed two-step method is demonstrably 
superior in resolving two incident angles that are quite close to 
each other. Fig.4(a) plots the correlation coefficient. The two 
angles are 0° and 25°. We note that the search region now 
range from —5° to 30°. Fig.4(b) illustrates the improvement 
achieved by using the proposed algorithm. 


Note also that the improvement in the resolution obtained 
by using the proposed approach is far superior to that realized 
when using the correlation or other traditional schemes. 
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Fig.2. For the case of simultaneous incidence of two waves 
coming from 0° and 3° (a) normalized magnitude; and (b) 
correlation coefficient. 
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Fig.3. For the case of simultaneous incidence of two waves 
coming from 0° and 3° (a) DOA results obtained by using the 
entire range; and (b) fine search DOA results (normalized). 
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Fig.4. For the case of simultaneous incidence of two waves 
coming from 0° and 25° (a) correlation coefficient; and(b) fine 
search DOA results(normalized) using the proposed algorithm. 
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